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Abstract. This work is motivated by the inverse conductivity 
problem of identifying an embedded object in an infinite slab. The 
novelty of our approach is that we use complex spherical waves 
rather than classical Calderon type functions. For Calderon type 
functions, they are growing exponentially on one side of a hyper- 
plane and decaying exponentially on the other side. Without extra 
modifications, they are inadequate for treating inverse problems in 
unbounded domains such as the infinite slab. The obvious reason 
for this is that Calderon type functions are not integrable on hy- 
perplanes. So they can not be used as measurements on infinite 
boundaries. 

For complex spherical waves used here, they blow up faster than 
any given positive polynomial order on the inner side of the unit 
sphere and decay to zero faster than any given negative polynomial 
order on the outer side of the unit sphere. We shall construct these 
special solutions for the conductivity equation in the unbounded 
domain by a Carleman estimate. Using complex spherical waves, 
we can treat the inverse problem of determining the object in the 
infinite slab like the problem in the bounded domain. Most impor- 
tantly, we can easily localize the boundary measurement, which is 
of great value in practice. On the other hand, since the probing 
fronts are spheres, it is possible to detect some concave parts of 
the object. 



1. Introduction 

Special solutions for elliptic equations or systems have played an im- 
portant role in inverse problems since the pioneering work of Calderon 
p. In 1987, Sylvester and Uhlmann |T3] introduced complex geo- 
metrical optics solutions to solve the inverse boundary value problem 
for the conductivity equation. Other type of special solutions, called 
oscillating-decaying solutions, were constructed for general elliptic sys- 
tems in PHI and ^2]- These oscillating-decaying solutions have been 
used in solving inverse problems, in particular in detecting inclusions 
and cavities |TU] . 

l 
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In developing the theory for inverse boundary value problems with 
partial measurements, approximate complex geometrical optics solu- 
tions concentrated near hyperplanes and near hemispheres for the Schr- 
odinger equation were given in [3] and jSj, respectively. In j^j, the con- 
struction was based on hyperbolic geometry and was applied in ^] 
to construct complex geometrical optics solutions for the Schrodinger 
equation where the real part of the phase function is a radial func- 
tion, i.e., its level surfaces are spheres. These solutions are called com- 
plex spherical waves in [J|. An important feature of these solutions 
is that they decay exponentially on one side of the sphere and grow 
exponentially on the other side. The hyperbolic geometry approach 
does not seem to work for the Laplacian with first order perturba- 
tions, such as the Schrodinger equation with magnetic potential and 
the isotropic elasticity. Recently, using Carleman estimates, complex 
spherical waves were constructed for the Schrodinger equation in [H], 
for the Schrodinger equation with magnetic potential in [2|, and for the 
isotropic elasticity in j!4j . The Carleman estimate is a more flexible 
tool in treating lower order perturbations. 

With these complex spherical waves at hand, one can study the in- 
verse problem of reconstructing unknown inclusions or cavities embed- 
ded in a body with known background medium. For the conductivity 
equation, this problem was documented in jl]. There are several re- 
sults, both theoretical and numerical, concerning the object identifica- 
tion problem by boundary measurements for the conductivity equation. 
We will not try to give a full account of these developments here. For 
detailed references, we refer to j3]. Here we want to point out that there 
are two striking advantages in using complex spherical waves for object 
identification problems. On one hand, one can avoid using a unique 
continuation procedure, more precisely, the Runge approximation, in 
the reconstruction. On the other hand, by the decaying property of 
complex spherical waves across the spherical front, the measurements 
can be easily localized. 

The construction of complex spherical waves and the inverse prob- 
lems studied in the literature cited above were restricted to bounded 
domains. However, there are many interesting inverse problems whose 
background domains are unbounded, for instance, identifying an object 
embedded in an infinite slab. In fact, our work here is motivated by 
this inverse problem. Ikehata in [S] studied the inverse conductivity 
problem in an infinite slab where the location of an inclusion is recon- 
structed by infinitely many boundary measurements. In he used 
Calderon type harmonic functions, i.e., e x i UJ + luJ± ) with uj G § ra_1 . These 
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functions are not integrable on hyperplanes. Therefore, they can not 
be the Dirichlet data of solutions with finite energy. To remedy this, he 
introduced Yarmukhamedov's Green function to construct a sequence 
of harmonic functions with finite energy that approximate the Calderon 
type function on a bounded part of the slab and are arbitrarily small 
on an unbounded part of the slab. 

In this paper we use complex spherical waves rather than Calderon 
type functions for the inverse problem in the slab. The most obvious 
advantage is that we do not need Yarmukhamedov's Green function to 
"localize" complex spherical waves since we can make these solutions 
decay faster than any given polynomial order on infinite hyperplanes. 
Therefore, we can treat the problem in the infinite slab like that in 
the bounded domain. Furthermore, we can handle the inhomogeneous 
background medium without requiring the Runge approximation prop- 
erty. Also, since our probing fronts are spheres, we are able to deter- 
mine some concave parts of the embedded object. 

This paper is organized as follows. In Section 2, we discuss a Carle- 
man estimate needed in our proof and its consequence. In Section 3, we 
construct complex spherical waves in the unbounded domain by means 
of the Carleman estimate. The investigation of the inverse problem in 
the infinite slab will be discussed in Section 4. 

2. Carleman estimate and related consequence 

As in j2], [H], and we construct complex spherical waves in the 
unbounded domain via a Carleman estimate. Let p G M. n , n > 2, be a 
fixed point. Let B £ (p) be the ball of radius e > centered at p. Without 
loss of generality, we take p = 0. We denote U = B £ (0) . From j2] and 
jH], we see that <p(x) = log|x| is a limiting Carleman weight for the 
semiclassical Laplacian —h 2 A in U. Namely, if a(x, £) = |£| 2 — \tp x \ 2 
and b(x, £) = 2<p' x ■ £, then 

{a, b}(x, £) = when a(x, £) = b(x, £) = 0. 

Here {a, b} = a'^ ■ b' x — a' x ■ b'^ is the Poisson bracket of a, b. With such 
(p, we have 

e -V/h _ M-l/h 

So e~ v l h becomes a polynomial weight. Carleman estimates with poly- 
nomial weights are well known in proving the strong unique continu- 
ation property for Schrodinger operators with singular potentials, see, 
for example, 0- Since we will work in L 2 -based spaces, it is enough to 
use a simpler L 2 Carleman estimate given for instance in [TT] . 
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Lemma 2.1. [TTl Lemma 2.1] There exists a C > such that for any 
v G C c °°(]R n \ {0}) and for any \ G {k + \ : k G N} with h < 1, we 
/iai>e i/iai 

/i 2 y |x|- 2 / ft |v| 2 |x|- n da; < C ^ |xr 2/,l+4 |/i 2 At;| 2 |xr n dx. (2.1) 

i/ere t/ie constant C is independent of v and h. 
Now we replace 1/h by 1/h — n — 5 such that 

h <E S := {(k + n + 5 + : G N}, 

where 5 > 0. Then implies that for t> G C^{U) and sufficiently 
small h £ S we have 

/i 2 y Ixr^ltflxr+^dx < C y \x\' 2 / h \h 2 Av\ 2 \x\ n+4+25 dx (2.2) 

with possibly different constant C. Let g G L™{U) and P ft * = |x|~ 1//l /i 2 (A- 
q)\x\ l l h . Then we get from ()2.2j) that 

^IMk^c/) < ^II^MIl| +2+4(C /) (2.3) 

for all v G C™(U) and small h <E S, where L 2 {U) is the weighted 
L 2 space with norm ||/||l2 (c/) = || (x) s f\\ L 2 (n) , where (x) = yjl + \x\ 2 . 
Combining (J2.3j) and the Hahn-Banach theorem, we prove the following 
existence theorem. 



Proposition 2.2. For any f G L_n_ s (U), there exists w G L 2 _n_ 2 _ 5 (U) 
such that 



P hW ■= \ x \ 1/h h 2 (A - q)(\x\~ l/h w) = f in U 

and 



h\H\ L 2 n {u) < C\\f\\ L 2 n {u) . 

Proof. Let us denote V = C™{U). For / G L 2 _n_ s (U), we define a 
linear functional £ on the linear subspace L = P£T> of L\ +2+5 {U) as 
follows: 

£(P*v) = (v,f) LHu) . 

Using (j2.3|) . we obtain 
\£(PhV)\ < \\v\\rt +s{u) \\f\\ L 2_^_ s(u) < ^||^IU| +2+5 ([/)ll/IU 2 _n_ 5 ([/)- 
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So by the Hahn-Banach theorem, £ can be extended to a linear func- 
tional on L\ +2 , 6 {U) with the same norm. Therefore, there exists a 

w G L 2 _n_ 2 _ s (U) such that for v E C c °°(fi) 

(P£V,W) L 2 {U) = (VJ) L 2 {U) 

and 

^HHIi^ (f/) < C||/|| L 2 n ([/) . 
The proof of proposition is now complete. D 

3. Complex spherical waves in unbounded domains 

In this section we will construct complex spherical waves for the 
conductivity equation in the unbounded domain Q = {x G MJ 1 : x n < 

c 

d} where d < 0. It is clear that Q C U(— B £ (0) ) for some e > 0. Let 
< j(x) E C 2 (U) and 7 = 1 in B R (0) C n SI for some R > 0. We look 
for a function v = Ix]' 1 ^!! with h E S satisfying 

Lv := V • (7VU) = in fi. (3.1) 

Using the Liouville transform w = ^/yu , we know that if w solves 

Pw := —Aw + qw = 

with <7 = Ay/^/^y then t> solves (|3.1j) . Here we note that q £ L£°(f2). 
We want to point out that by checking the proofs in jj] , complex spher- 
ical waves constructed there can be extended to the unbounded domain 
Q. In this work we take a different route using the Carleman estimate 
which can be adapted to other equations more easily. 

As in j2j and [Hj, we want to construct solutions to Pw = in Q 
which have the form 

w = e- p{x)/h (a + r), 
where p = ip + iip and <p, ip are given by 

<p(x) = log \x\, 
ip(x) = dist§n-i(|||,ei). 

Our first goal is to see what a will look like. To this end, we change 
coordinates and write ') where x' = r8 with r > and 

6 E § n " 2 . We take z — X\ + ir to be a complex variable, and let 
\1/ : x i— > (z, 9) be the corresponding change of coordinates in Q. If / is 
a function in Q we write / = / ° Thus we can see that 

p = log z } 

(V p)~ = - (e\ + ie r ) with e r — (0,9), 
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and 

z\z — z) 

Also, Vp ■ V becomes (2/z)dz in the new coordinates. 
Writing h 2 P = (—ihV) 2 + /i 2 g, we have 

e p/h h 2 Pe- p/h = {-ihV + zVp) 2 + h 2 q 

= -(Vp) 2 + h(2Vp • V + Ap) + h 2 P 

= h(2Vp ■ V + Ap) + h 2 P. (3.2) 

Here, <p and ip were chosen so that (Vp) 2 = 0. To get P(e~ p ^ h (a + r)) = 
0, we choose a and r to satisfy in Q 

(Vp ■ V + ^Ap)a = (3.3) 

and 

e p/h h 2 Pe- p/h r = -h 2 Pa. (3.4) 

The first equation ()3.3|) is a transport equation for a. Writing a = ao$, 
this reduces to 

n — 2 

( * " 2<73iy )a = »*<«>■ 

It is easy to see that the general solution is given by 

a = (z - Z) 2 g 

for any g(z,9) satisfying d z g = 0. 

Taking g = 1 and going back to the x coordinates, we obtain 

2 — n 

a(x) = (2i\x'\) — , xeVL. (3.5) 
In (13. 4|) . the right hand side becomes 

o > 2 — n o / 1 , v. >. — 2 — n ,, 2 — n \ 

-/i 2 Pa = (2i)T- /i 2 ^-(n - 2)(4-n)|x'| _ 3- - — J. 

This decays in x' but is constant in xi, so Pa G £ 2 (^) for s < — n/2. 
We now rewrite (|3.4|) as 

Since e~ l ^^ h h 2 Pa £ L 2 (il) with s = —n/2 — 5 for some 5 > 0, we deduce 
from Proposition 12.21 that there exists r solving (J3.4)) and satisfying 

Mb„ „ M)<Ch. (3.6) 



Therefore, we have constructed the special solution w of — Aw + qw = 
in Q having the form 

w = e -(<fH*i>)/h( a + r ) = \ x \- x l h e-^l h (a + r) 



COMPLEX SPHERICAL WAVES IN UNBOUNDED DOMAINS 7 

with a given by (j3.5j) and r satisfying (jH.fijl . In particular, we can 
see that u> G L 2 (f2) for all sufficiently small ft,. Furthermore, since 
Aw = qw G £ 2 (fi), we can use a cut-off technique and the elliptic 
estimate to show that 

w G H 2 (n), 

where Cl = {x G R n : x n < d} for any d < d. In other words, the trace 
of w on any hyperplane Hf, := {x„ = b} with 6 < d is well-defined and 

d Xn w\ Hb e H 1 ' 2 . 

This property is an important difference between complex spherical 
waves and Calderon type solutions. 

In order to apply complex spherical waves to the inverse problem, 
we need to estimate the H 1 norm of r on any fixed open subset of Q. 
Let fli and Q2 be two open bounded subsets of Q such that Q\ C ^2 
and ^2 C fi. From (|3.2j) we see that 

e p/h h 2 Pe- p/h = -h 2 A + h(2Vp ■ V + Ap) + h 2 q. 

Therefore, e p/h h 2 Pe~ p/h is elliptic as a semiclassical operator. Further- 
more, we have 

\\e p / h h 2 Pe- p / h r\\ L 2 {n2) = ||/i 2 Pa|| L2(f , 2) < Ch 2 

and 

lklU 2 (n 2 ) < Ch. 

Hence, by the semiclassical version of elliptic estimate (HJ Lemma 2.6], 
we obtain that 

Wd a r\\ L * m < Ch, 

\a\<2 

in particular, 

IMU»(no + ^||Vr|| L 2 (ni ) < Ch. (3.7) 

Finally, let p be any point such that dist(p, Q) > 0, then complex 
spherical waves for (J3.1)) in Q are given by 

v = v (x, h) = 7 ' 1/2 |x - V \- l l h e~^ th {a + r) 

2 — n 

for h G S small, where if) = dist§n-i (jfz^f 5 e i)' a ( x ) = \Zi\x' — p'\)^~ 
and r satisfies (j3.6|) and (j3.7|) . Note that for all x G Q we have x' — p' = 
(x 2 - p 2 , ■ • • , x n - p n ) ^ since p n > x n . 
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4. Inverse problem in a slab 

In this section we shall apply the special solutions we constructed 
above to the inverse problem in an infinite slab. We study the recon- 
struction of an embedded cavity here. The same method works for the 
inclusion case and so does the method in [3] or We leave this 

generalization to the interested reader. Let Q = {x G M. n : d\ < x n < 
d 2 , di < d 2 }. Assume that D is a domain with C 2 boundary in Q so 
that D C VL and Q \ D is connected. Let u(x) be the unique solution 
of finite energy to 

L 7 w = in Q\D, 

7 §h = on 3D, (4.1) 
u = f G H 1 ' 2 on dQ, 

where the conductivity parameter j(x) satisfies the assumptions de- 
scribed in the previous section. The well-posedness of (14.1)1 can be 
proved by the standard Lax-Milgram theorem with a Poincare-type in- 
equality in the infinite slab: there exists C > independent of u such 
that 

I \u\ 2 dx <C I \Vu\dx (4.2) 
Jn\D Jn\D 

for all u G H l (VL \ D) with u = on x n . = d\ and di. One possible way 
to prove (|4.2|) is to divide Q = Z U (Q \ Z), where Z is a cylindrical 
domain of the form Z = B x (d%, d 2 ), B is a ball in M n_1 , and D C Z. 
Then (14.2)1 is a consequence of 

tfdx < C I \Vv\dx (4.3) 
z Jz 

for all v G H l (Z) with v = on <9Z fl {x n = di, d 2 } and 

/ \w\ 2 dx<C [ \Vw\dx (4.4) 

for all w G H\Q\Z) with w = on d{Q\Z)n{x n = d 1 ,d 2 }. Now (Oil 
can be proved by a contradiction argument and (|4.4)l can be established 
by the usual integration technique. For brevity, we left the details to 
the reader. 

The inverse problem here is to identify D from the Dirichlet-to- 
Neumann map A D : H^ 2 (dQ) i-> H-^ 2 (dQ) 

Ad : / >-> 7n~ an- 
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In the weak formulation, the Dirichlet-to-Neumann map is defined by 



(A D f,g) = / 7Vw Vvdx 
Jn\D 



where v G H l (Q \ D) with g = v\q^. We are interested in the recon- 
struction problem in this work. 

Let Q = {x G 1R 3 : x n < d} with g?2 < d and t > 0. Pick any point 
p so that dist(p, ft) > 0. We now denote complex spherical waves of 
LyV = in Cl by 



^(a;, /i) = ^A 7 -l/2| ;r _ p |-l//, e -^A (a+r) = 7 -l/2 ( ^^ ) l//^ e -^^( a+r ). 

\x — p\ 



Let us define the energy gap functional 



E t (h) := ((A -A D )v t (x,h),v t (x,h)) = / (A -A D )f t (x, h)-v t {x, h) ds, 

Jan 



where A is the Dirichlet-to-Neumann map for L 7 with D = 0. Note 
that E t (h) is well-defined for any £ > even when dfl is unbounded. 
To understand how we reconstruct D, we first observe that E t (h) can 
be estimated by 



1 

C 



\Vv t (x,h)\ 2 dx<E t (h)<C / (\Vv t (x,h)\ 2 + \v t (x,h)\ 2 )dx 

D J D 



(4.5) 

for some constant C > 0. 



Proof of ([4.5)1 . Let u be the solution of (j4.1|) with / = t> t . We observe 
that 

/ 7VM ■ V(w - u t )da; = 0. (4.6) 
Jn\D 
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Using the definition of the Dirichlet-to-Neumann map and (|4.fij) . we 
can compute 

((A -A D )v t ,v t ) 

dv t f du 

7— — v t ds — I ry—Vtds 



on du J dn dv 



~jVv t ■ Vv t dx — I •y'Vu ■ Vv t dx 

U\D 



"fVv t ■ Vv t dx — / 7VU • Vv t dx + / 7VM • V(« — v t )dx 
n\D Jn\D Ju\D 



+ / i\Vv t \ 2 dx 
Jd 

= / 'j\Vv t \ 2 dx+ / 7V(u - v t ) ■ V{u - v t )dx. (4.7) 
Jd Jn\D 

The formula ()4.7j) immediately implies the first inequality of (|4.5jl . 

To obtain the second inequality of (|4.5jl . we need to estimate the last 
term in (|4.7jl . Denote w = u — v t . By the construction, we have that 

L 7 u> = in tt \ D, 

7f = -7fe on 0D, 
w = on <9fi. 

It follows from the elliptic regularity theorem that 

\H\m(n\D) < C\\i^\\ H -i/2 {dD) . (4.8) 
On the other hand, we know that L 7 vt = in D and therefore, 

h^\\H-i/*(dD) < C\\v t \\ H i {D) . (4.9) 

Putting together (|4~£|) and leads to 

IMItfHnYD) - C W V t\\H^(D) 
and the second estimate in ([4.5)1 . D 

Now we will use (|4.5jl to study the behaviors of E t (h) as ft — > and 
h E S for different t's. 

Theorem 4.1. Let dist(7J>,p) =: c?o > 0. For t > and sufficiently 
small h £ S, we have that 

(i) . 7/ d > t, then E t (h) < Ca 1/h for some a < 1; 

(ii) . If d Q < t, then E t (h) > C(3 l/h for some (3 > 1; 
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(iii). IfDHBtip) = {y}, then C- l h n ~ 2 < E t {h) < Cbr 1 for some 
C>0. 

Proof, (i). If dist = d > t then 



\x — p\ 



<1 VxeL>. 



Combining the second inequality of ()4.5|) and the behavior of v t (x, h), 
we immediately prove this statement. 

(ii). We first pick a small ball B s CC B t (p) n D. Using (JH2D, the 
leading term of Vv t (x, h) is 

1 f t \l/h„-ty/h„-l/2 



h \x — p\ 



7 _i/i V(log |x -p| +iif))a. (4.10) 



We notice that 7 1//2 V(log \x — p\ + iip)a 7^ for all x E Bs- So this 
statement follows from the first inequality of (|4.5|) and the fact 



, 1 , > 1 V x G £ a . 

|x — p\ 

(iii). Pick a small cone with vertex at y, say T, so that there exists an 
e > satisfying 

T n {0 < \x - y\ < e} C D. 

It is not restrictive to take y = and p — (0, • • • , 0, t) for t > 0. Now 
we observe that if z e T and |z — y| = \z\ = s < e then 

\z — p\ < s + t, 

that is 

t t 



> 



\z — p\ s + t 

Hence, from the first inequality of (|4.5j) and (|4.1(J|) . for h 1 we have 
that 



> 

~ h 2 ^2-kh 

> Ch n ' 2 . 



On the other hand, we can choose a cone T with vertex at p such 
that D C T fl {t < \x — p\ < t + r]} for 77 > 0. Thus, by the second 
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inequality of (|4.5jl and (|4.1()jl . we can estimate 




< ch-\ 



□ 



The distinct behavior of E t (h) immediately allows us to detect the 
boundary of the cavity. However, Theorem 14.11 is unpractical since 
we need to take the measurements on the whole unbounded bound- 
ary. Fortunately, taking advantage of the decaying property of com- 
plex spherical waves, we are able to localize the measurement near the 
part where complex spherical waves are not decaying. To be precise, 



where 5 > is sufficiently small. Now we are going to use the mea- 
surement fs t t(x,h) = <f>s,tVt(x,h)\dn^ Clearly, the measurement fs,t is 
localized on B t +s{p) H dQ. In fact, if 5 is sufficiently small, fs,t is lo- 
calized on only one part of dQ, depending on whether p lies above or 
below Q. Let us define 



Proof. The main idea is to prove that the error caused by the remain- 
ing part of the measurement gs,t '■= (1 — 4>s,t)vt(x, h)\gn is as small as 
any given polynomial order. Let ws,t(x, h) be the unique solution of 



with boundary value gs,t- We now want to compare wg )t with (1 — <fi$j)v t . 
To this end, we first observe that 



I L 7 ((l - (f>6,t)v t - w s , t ) = £ 7 ((1 - 4>6,t)v t ) in Q, 
\ (1 - 4>s,t)v t - ws,t = on d£l. 

Notice that 

L 7 ((l - (f) S> t)v t ) = (1 - (p 5 ,t)L^v t + [L 7 , (1 - <f)s,t)]v t = [L 7 , (1 - <f>s,t)W 
So we see that 

supp (L y ((l - (p s ,t)v t )) C H n {t + 5/2 < \x - p\ < t + 5} 



let 4> Sit (x) e C °°(M n ) satisfy 
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and therefore 

||L 7 ((1 - <f>s,t)vt)\\i?(p) < Ca\ /h 
for some < a\ < 1. Consequently, we have that 

||(1 - <j>s,t)v t - w s , t \\m(a) < Ca{ /h , 

in particular, 

||(1 - Mv t - w s ,t\\m { D) < Ca\ /h . (4.11) 

Using (|4.5jl for ((A — Ao)g5,ti gs,t) with v t being replaced by wg jt , we 
get from (|4.11|) and decaying property of v t that 

((A - A D )g s , t , gs,t) < Cal /h 

for some < a 2 < 1. 

Now we first consider (i) of Theorem 14.11 for E^ t {h). We shall use a 
trick given in [I]. In view of the definition of the energy gap functional 
and the proof of the first inequality of ([4.5)1 . we get that 

o < <(A - A D )(CA t ± C%,t), Us,t ± C l gs,t) 

for any ( > 0, which leads to 

|((A - A D )fs,t,96,t) + <(Ao - A D )# M , A*) I 
< C 2 ((Ao - Ad) At, fs,t) + C 2 ((A - A D )g s , u 9s,t). (4-12) 

It now follows from v t (x, h)\gn = fs,t + 9s,t and (|4.12j) with £ = l/\/2 
that 

-((A - A D )fs,t,fs,t) 

< ((A - Aj?)&,t, # 5)t ) + ((A - A D )% u t ) 

< C^ //l + ((A -A D K,^). (4.13) 

So from (i) of Theorem 14. 1[ the same statement holds for Eg )t (h). Like- 
wise, the second inequality of (iii) in Theorem 14.11 also holds. 

Next we consider (ii) and the first inequality of (iii) in Theorem 14.11 
for Es t t- Choosing £ = 1 in ()4.12|) we get that 

^((A -A D )v t ,v t ) 

< ((A - A D )g s ,t, gs,t) + ((Ao - A D )f s ^ f 6jk ) 

< Ca\ /h + ((A -A D ) fs, t J 5>t ). (4.14) 

Therefore, (ii) of Theorem 14.11 and ()4.14|) implies that the same fact 
holds for Es t t- Likewise, the first inequality of (iii) in Theorem 14. 11 holds 
for Es t t- The proof is now complete. LJ 



14 



SALO AND WANG 



Remark 4.3. With the help of Theorem \4- 6 <\ when parts of 3D are 
near the boundary dVL, we could detect some partial information of dD 
from only a few measurements taken from a very small region on one 
side ofdQ. 

To end the presentation, we provide an algorithm of the method. 

Step 1. Pick a point p Q and near dQ. Construct complex spherical 
waves v t (x, h) for h G S. 

Step 2. Draw two balls B t (p) and B t +s(p). Set the Dirichlet data 
fs,t — <i>8,tVt\dn- Measure the Neumann data Apfsj over the region 

B t+s ( P )ndQ. 

Step 3. Calculate Eg t t(h) = ((A — Arj)fs,t,fs,t)- If Eg,t(h) tends to 
zero as h — > 0, then the probing front {|se — p\ — t} does not intersect 
the inclusion. Increase t and compute E^ t (K) again. 

Step 4. If Es t t(h) increases to oo as h — > 0, then the front {\x— p\ = t} 
intersects the inclusion. Decrease t to make more accurate estimate of 
dD. 

Step 5. Choose a different p and repeat the procedures Step 1-4. 
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